On the Relation between Optical Conductivity and Quasiparticle Dynamics: Boson 

Structures 
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An extended Drude form is often used to analyze optical data in terms of an optical scattering rate 
and renormalized mass corresponding, respectively, to the real and imaginary part of the memory 
function. We study the relationship between memory function and quasiparticle self energy for an 
isotropic system. We emphasize particularly boson signatures. We find it useful to introduce a new 
auxiliary model scattering rate and its Kramers-Kronig transform determined solely from optics 
which are much closer to the self energy than is the memory function itself in the normal state. In 
the superconducting state the simplification fails because the quasiparticle density of states acquires 
an essential energy dependence. 

PACS numbers: 74.20.Mn 74.25.Gz 74.72.-h 



I. INTRODUCTION 



Opticalii^ and angular resolved photo emission 
(ARPES) dataiSi^ have given a wealth of information on 
quasiparticle dynamics in the cuprates both in their nor- 
mal and superconducting state. The methods are com- 
plementary but the exact quantitative relationship be- 
tween the two is complex. One important difference is 
that ARPES gives direct information on angular varia- 
tions around the Fermi surface while optics involves an 
average over all the quasiparticles participating in the 
absorption. Even if these anisotropics are not accounted 
for (isotropic system) there remain additional differences 
which have their fundamental origin in the fact that 
ARPES measures directly the quasiparticle spectral den- 
sity A(k, Lo) at fixed quasiparticle momentum k as a func- 
tion of energy ui while optics involves the current-current 
correlation function which depends on the product of two 
spectral densitiesSii at the same momentum k but with 
frequencies displaced by the photon energy. 

In this paper we focus on similarities as well as on 
essential differences between the information that is de- 
rived from these two different probes and particularly on 
how they are to be compared. We will be interested in 
both, the normal and the superconducting state at zero 
and at finite temperature. In the literature on optics it is 
the temperature and frequency dependence of the opti- 
cal scattering rate T~p(T, oj) which has been particularly 
emphasized)^ More recently the imaginary parli^ of the 
memory functior&iS related to T^p^{T,uj) by Kramers- 
Kronig (KK) transform has also been used to compare 
directly with the energy dependence of the real part of 
the quasiparticle self energy determined by ARPES>ii*i^ 
Of particular interest in such a comparison is the rela- 
tionship between the position in frequency of peak struc- 
tures seen in these quantities and how they reflect cor- 
responding peaks in the electron-boson spectral density 



/^x(w)iiii^ii^ For phonons J^x(tj) is the well known 
electron-phonon spectral density while in the cuprates 
exchange of spin fluctuation32iiSiiLi& is more appropriate 
although as yet there exists no consensus— on this issue. 
Here we will consider several models for I'^xi^) within an 
extended Eliashberg formalism. In the superconducting 
state provision is made for d-wave symmetry of the gap. 
To keep things simple the same form of /^x(tj) is used in 
both gap and renormalization channels but with different 
magnitudes with the ratio between the two equal to g. 

In Sec. |n]we consider only the normal state and cou- 
pling of quasiparticles to a single {uje) Einstein mode&SS 
We also employ a simplified approximate scattering time 
formulation of the optical conductivity formula in the 
normal state. This allows some analytic results to be es- 
tablished and provides the motivation for introducing a 
new model scattering rate and its KK-transform deter- 
mined solely from optical data but which is very close 
to the self energy itself. While the real part of the self 
energy has a logarithmic singularity at loe the imaginary 
part of the memory function does not. Instead, it shows 
a peak at w = v^w^- On the other hand, our model 
quantities reproduce well the self energy at T = and 
only very small differences arise at finite T. In Sec. IIIII 
we provide the more general formalism needed to treat 
the conductivity accurately as well as to take care of ex- 
tended Px{^) spectra and the superconducting state. 
Results in these cases are presented in Sec. IIVI For the 
superconducting state additional complications arise be- 
cause of the essential energy dependence of the self con- 
sistent quasiparticle density of states (DOS). In Sec. IVl 
we present our numerical results for our new model scat- 
tering rate based on the complete formula for the op- 
tical conductivity and solutions of the Eliashberg equa- 
tions and compare these with the imaginary part of the 
quasiparticle self energy. We confirm the previous ex- 
pectation, based on simplifying assumptions, that in the 
normal state these two quantities are close in magnitude 
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as well as frequency dependence. In the superconduct- 
ing state they are not. Finally, conclusions are found in 
Sec. ED 



II. NORMAL STATE AND COUPLING TO A 
SINGLE BOSON 



A generalized Drude form with frequency dependent 
optical scattering rate t~^{ijj) and optical effective mass 
ratio [m* (w) / m]op has been used extensively to analyze 
optical data in a one component approach. The opti- 
cal conductivity (t(w) is written in terms of the mem- 
ory function^ M{u)) = t^^{lo) — iu}Xop{uj) with Xop{cj) = 
{[m*{uj)/m\op — 1} the optical mass renormalization pa- 
rameter. We havei^ 



cr(w) — cri(tj) + 10-2 (w) 



P 



47r iTop^{Lu) +uj[l + Xopiuj)] 

(1) 

and Tgpluj) and iujXop{uj) are related by Kramers Kro- 
ning, so that when the optical scattering rate is known, 
the optical mass renormalization can be obtained from 
the KK-transform. 

The conductivity is a two particle property given by 
the current-current correlation function. It is related in 
a complicated way to the quasiparticle self energy (a one 

In the normal state-- at zero 



particle property) S(a;) 
temperature (T = 0) 



A-KU) J , 
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(2) 



where t"^^^ is the impurity scattering rate and S(i') 
accounts for the inelastic scattering. In terms of 
the electron-phonon or electron-spin fluctuation spectral 
density a'^F{uj) or Pxi^) respectively2i2^ [call it F{uj)] 
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FIG. 1; Comparison of the real part of the quasiparticle self 
energy Ei (uj) as a function of u) with the corresponding optical 
quantity uj{[m'' (uj) /m]op — 1} = ujXopii^)- Coupling is to a 
single Einstein oscillator with uje = 40meV and A = Atr = 
1 meV in Eqs. |0J and ||5J. There is a logarithmic singularity 
in Ei(ti;) at lj = uje while uj\o'p{lo) only has a small peak 
at ci; = \/2uiE and a logarithmic singularity in the slope at 

LO — UJe- 



with the Greens function G(k, w -I- iO+) = [tj — Ek — 
+ jO+)]~^, where ek is the quasiparticle dispersion. 
A useful! approximate expression for T~p[uj) was ob- 
tained directly in second order perturbation theory by 
P.B. AUemSa It is (for T = 0) 
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dzF{z){uj- z) + Ti„lp. 



(6) 



Details in the validity of this form, based on numerical 
evaluation of Eqs. (0) and (O are found in Refs. and 
m For > 
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dnFifi). 



(3) 



For coupling to a single Einstein mode F(uj) = A5{lo 
cj^;) we obtain the simple, well known formula 



Y.{uj) = Aln 



UJe 



UJe 



iiTA9{\uj\ — UJe) 



(4) 



with d{x) the Heavyside 0-function. It is clear from Q 
that the imaginary part has a finite jump at uje as well 
as a logarithmic singularity in its real part which can be 
used to identify boson structure in both cases. The self 
energy can in principle be determined from ARPES data 
which measures the quasiparticle spectral function 



A{k,uj) 



1 



IniG(k,w-f iO+), 



(5) 



where the subscript tr means transport, although for sim- 
plicity, here we will not make a distinction between quasi- 
particle and transport spectral density; in general they 
are different. (We have set the impurity term equal to 
zero for simplicity.) Note that this formula is less favor- 
able for the identification of boson structures than is the 
imaginary part of the self energy. While it is zero for 
UJ < UJe and finite for uj > u>e it has no jump at uje- 
Instead, it smoothly increases from its zero value which 
makes it harder to identify the exact position oi uje in 
T~p (w). AppUcation of the KK-transform to Eq. 10 gives 
immediately 



ujXop{uj) ~ —2A In 



UJe 



UJE — UJ 



— in 



(8) 



as the optical mass enhancement. This is to be compared 
with the real part of the quasiparticle self energy, Eq. Q). 
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In terms of the usual mass enhancement parameter A 
[defined as two times the first inverse moment of F{uj)], 
A — Xlue/'2, and the limit of Si(a;) as a; ^ equals 
—a; A. Thus the slope of the self energy gives the value of 
A as does also cuXop ~ Xuj. This is seen in Fig. ^ where 
we compare Si (a;) (solid line) and LuXop{uj) vs uj (dashed 
line) for parameters ^ = 1 meV and w^; = 40 meV. While 
these functions agree in the w — > limit they deviate from 
each other at finite frequencies. In particular, at w = uje 
the quasiparticle self energy exhibits a logarithmic sin- 
gularity while the optical mass renormalization coXop^Lo) 
has only a logarithmic singularity in its slope. It can, 
furthermore, be shown that uXopito) has its maximum at 
UJ = V2uje- Thus, while the boson structure in T,{uj) ap- 
pears at the boson energy and is singular, the maximum 
in Aop(w) is instead displaced to \/2uje and is associated 
with a rather broad peak in comparison. We conclude 
from these considerations that the optical mass enhance- 
ment parameter Xopito) can be used to identify boson 
structure in the imaginary part of the memory function 
as done recently by Hwang et al? but its signature is 
much weaker than is the case in the quasiparticle self en- 
ergy and w_B is displaced by a factor of v^, a fact that 
does not appear to have been appreciated before, but is 
important to realize. 

In preparation for what will come later when we con- 
sider the superconducting state, we note that the self 
energy Q and, consequently, the corresponding opti- 
cal quantity JSJ would also applji2i*22iS^ in an impu- 
rity model with constant quasiparticle scattering rate 
T~p = A, but in which the normal state quasiparticle 
density of states, N{e), has a gap of size we at the Fermi 
surface, i.e.: N{e) = Ni^ ior e < Q and e > we but is zero 
for Q < e < loe- This shows immediately that structure 
in N{e) can have an effect on the self energy which is, in 
some cases, indistinguishable from boson structure. We 
return to this important point later. 

How is this picture changed when we consider finite 
temperatures? In this case the formulas determining 
the conductivity as well as the quasiparticle self energy 
are more complex. Nevertheless, a simple picture still 
emerges. The normal state optical conductivity is now 
given bji2& 
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dv tanh 
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^ S-^{T,uj,v) 



dv 



tanh ( ^ \ - tanh 
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xS-\T,uj,v)^, 



(9) 



where 



S{T, uj,y)^uj + E*(r, v + Lo)- I](T, y) - rr^^, (10) 



and the self energy ii 
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where tpi^) is the digamma function and S*(r, w) is the 
complex conjugate of I](T, tj). For a general F{z) 



T^}{T,uj) = -2S2(r,w) 



TT ldzF{z) 

' UJ + z 



2coth ^ 



— tanh 
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(12) 



Shulga et alm^ were able to show that to a good approx- 
imation 



Top^{T,i^) 



- (dzF(z) 
^ J 

— (cj + z) coth 
~\~{uj — z) coth 



2uj coth I — - 
2T 



UJ + z 



2T 

UJ — z 

2T 



(13) 



At zero temperature Eqs. (|12|l and (|13|) reduce to the 
imaginary part of Eq. ^ and to Eq. ©, respectively. 
KK-transforms can be used to get from Eqs. H12|l and l|13|l 
the real part of S(T, uo) at finite T as well as ujXop{T, uj). 
Results are shown in Fig.|21for coupling to a single Ein- 
stein mode as in Fig. ^ Three curves are show in each 
frame. The dotted line applies to the quasiparticle self 
energy, the dashed line to the corresponding optical quan- 
tity and the solid line to a model yet to be defined. The 
top frames give real parts at T = Q.Iuje (left) and at 
T = 0.25uje (right). The corresponding scattering rates 
on which these are based, are shown in the two bottom 
frames. First we consider the top frames. In all cases 
we see that temperature smooths out the correspond- 
ing boson structures, but even for T = 0.25uje they re- 
main easily identifiable, although, the difference between 
quasiparticle and optical quantities is no longer as pro- 
nounced. 

The solid curve in Fig. [5] is based on the following ob- 
servation. At T = 0, d[u}T~p^{uj)]/duj given by Eq. ® 
equals exactly the quasiparticle scattering rate given as 
twice the imaginary part of the quasiparticle self energy 
in Eq. Thus, in this particular limit one can get the 
quasiparticle scattering rate directly from optics by tak- 
ing the first derivative of ujT^p{uj). Of course, we have 
assumed that anisotropics in momentum space can be 
neglected. In general, ARPES will give information on 
the k-dependence of the quasiparticle self energy while 
optics is always an average. Neglecting such complica- 
tions, it is interesting to introduce a model scattering 
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FIG. 2: Top frames give a comparison of the real part of the 
quasiparticle self energy 2Ei(ij) (dotted line) with its optical 
(dashed line) and model (solid line) counterparts as a function 
of LO. The right hand frame is for temperature T = 0.25aj_B 
and the left hand frame is for T = O.Iljb. The bottom frames 
gives corresponding quasiparticle scattering rates on which 
the results in the two top frames are based. The parameters 
are A — 1 meV and cjb = 40 meV. 



rate t^^^^j = d[ujT^p^ (lu)] / duj defined from optics alone 
and consider its relationship to the quasiparticle self en- 
ergy at finite temperatures. The formula for t~],^^i{T,u!) 
equivalent to Eqs. ((T^ and ((T^ is 



model 



{T,uj) = TT / dzF{z) 

4(^) 



2coth 



+ 



(2r) 

-4(^) 

sinh^ sinh^ {^) 



1 



.(14) 



In all cases considered in Fig. [21 the model quantity 
(solid lines) deviates from its quasiparticle counter part 
(dotted lines) only very slightly while it deviates much 
more from the optical case (dashed lines). It is clear, 
therefore, that if one wishes to compare optical and quasi- 
particle quantities directly, it is better to use the KK- 
transform of t~1^^j{u;) than that for t~J-{lj) itself. The 
peaks in the corresponding model mass enhancement will 
now correspond to the Einstein frequency loe and it will 
be sharper than in the optical quantity itself which, as 



we have seen, is also shifted by \/2- 

Note also, that on comparing top and bottom frame 
it is the function t^^ (cj) which is least useful in localiz- 
ing boson structures as discussed in relation to Eq. 
All the other functions have a sharper signature of the 
Einstein frequency uje- 

We next turn to a discussion of how the use of an ex- 
tended spectrum for F{lu) modifies our simple results and 
how the superconducting condensation further modifies 
boson structures in quasiparticle and optical quantities. 



III. EXTENDED BOSON SPECTRUM AND 
SUPERCONDUCTING STATE 



In general, coupling of the electrons to a boson spec- 
trum such as phonons or spin fluctuations is not re- 
stricted to a single mode. The corresponding a'^F{uj) 
(phonons) or I'^xi^) (spin fluctuations) can extend to 
100 meV or so and even up to 400 meV (of the order of J 
in the t ~ J model^) , respectively. Here we wish to con- 
sider such extended spectra and also consider the super- 
conducting state assuming d-wave gap symmetry since 
we have the cuprates in mind. The equation for a{uj)^ 



2 1°° \ 
air, v)^^l.ljduj tanh (^^) [J(c., v) ~ J{-u, v)] \ 

(15a) 

where {■ ■ ■)e denotes the averaging over angle 9 and the 
function J{lo, v) is given by 



^^^^''^^ = m^F^ ^ [l^iV(^;0) 

xN{uj + v;e)- P{lo; 0)P{uj + v; 9)] 
1 

^ E*{uj]9)- E{Lu + v-e) 
X [1 + N*{uj;e)N{uj + v;e) 
+P*{uo-e)P{Lu + v-9)], (15b) 



with 



E{uj; 0) = + iO+) ~A^{uj + iO+; 9), (15c) 



and 



N{uj;0) 



oj{lj + iO+) 
E{LO-e) ' 



p{u-e) 



A(cj + iO+; 



(15d) 

Here E* , N*, and P* are the complex conjugates of E, 
N, and P, respectively. The two fundamental functions 
Lu^uj) and A(w) are related the renormalization and gap 
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functions, respectively^ They are for the gap channel 
A(i/ + iO+;6') 

OC 

= ^Tg J2 cos(20) [A( 



0j^{iuJrn) + A'^iiuJm; 0') , 





+ ing I dz cos(26')/^x(z) [n{z) + f{z - i^)] x 

A(i/-z + iO+;6i')cos(26l') \ 
{ly- z + iO+) - A2(i/ - z + iO+; 61') /, 

and in the renormalization channel 

+ i6) 

oo 

— ly + inT [X{iy — iuim) ~ + ium)] 

m=0 



(16a) 




a)2(iw„) + A2(icj„; 6') I 
+ iTr I dzl'^x{z)[n'{z) + f{z - i^)] 

Lj{i'-z + iO+) 
{u-z + - A2(z/ -z + iO+; 61') , 



Here 

The gap is given by 
A(j/ + iO^ 



A{iy + i0+;9' 
oj{iy + iO+) 



(16b) 



(16c) 



(16d) 



or, if the renormalization function Z{iy) is introduced in 
the usual way as a)(i^ + iO^) = iyZ{v) then 



A{iy + iO+; 



(16e) 



These equations are a minimum set and go beyond a BCS 
approach. They include inelastic scattering known to be 
strong in the cuprate superconductors. In Eq. p6a|l g 
modifies the electron-boson spectral density I'^xi^) from 
its value in the renormalized channel. In general the 
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FIG. 3: Top frame: the three models for the electron-boson 
spectral density I'^xi'^) used in Fig.|l| The MMP model with 
ijsF = 41 meV, Xmaip ~ 0.554 augmented by a (5-function 
at uiE ~ 41 meV with Xs = 0.946 for a total of 1.5 (dashed 
curve). The model spectrum with A = 2.01 obtained from a 
fit to the optical conductivity as described in Ref. (dot- 
ted curve). Coupling to an optical resonance at 41 meV is in- 
cluded as a peak. A pure MMP spectrum with uisf = 41 meV 
and A — 1.5 (solid curve). Bottom frame: the superconduct- 
ing density of states iV(a;)/iV(0) at T = lOK for the three 
spectra shown in the top frame. 



shape of I'^xi^) could also be different but we have not 
included this possible complication here. 

In what follows, we will use for I'^xi'^) ^ form suggested 
by Millis et al<^ (MMP) given by 



lo/losf 



(17) 



SF 



with ws F a spin fluctuation frequency which can be fitted 
to optical data and to emphasize structure we also will 
add in one case a (5-function at some specific frequency. 
Such an MMP form for losf — 41 meV is shown in the 
top frame of Fig. |31 (sohd fine). It has a A = 1.5. Also 
shown is the same MMP form now with an added i5- 
function peak at loe = 41 meV (dashed line). A is again 
equal to 1.5 but now with only 0.554 in the MMP form 
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and the rest in the 5-function. In the superconducting 
state the MMP form is modified because of the growth of 
an optical resonanceiS'ii'i^ at 41 meV which is the same 
frequency at which a spin one resonance is also observed 
in spin polarized inelastic neutron scattering2^ for opti- 
mally doped YBCO. (In other cuprates the position of 
the observed resonance varies somewhat j2^»24) This opti- 
cal resonance grows as the temperature is lowered with 
the area under the resonance scaling as the superfluid 
densitjii^ to a reasonable approximation. Such a spec- 
trum has been derived from a fit to T = 10 K optical data 
on a YBa2Cu3 6.95 (YBCOg.gs) twinned single crystal. 
This spectrum was also used later on to fit optical data 
for untwinned single crystals reported by C.C. Homes et 
al^^^ and also to calculate the microwave conductiv- 
ity of the YBCOg.gg single crystals. ^^'^^ This spectrum 
is shown as the dotted curve in the top frame of Fig. 
It has a A = 2.01. We will take this I'^xi^) spectrum 
as representative of the oxides. Note that coupling to 
a collective mode at (tt, tt) is also seen in the ARPES 
data of Campuzano et ali^ Also note that in the above 
formulation we have neglected possible momentum space 
anisotropics. Recent ARPES data by Kaminski et a/.'^° 
justify this assumption. 

The bottom frame of Fig. El shows the quasiparticle 
density of states N{uj) /N{0) in the superconducting state 
for T = lOK for the three spectra presented in the top 
frame. This would be the classic way to see boson struc- 
ture in conventional superconductors through tunneling 
which, so far, has been less successful in the high Tc 
oxides. This is the reason why optics and ARPES be- 
came such important tools in studying the quasiparti- 
cle properties. Note that the solid curve for N{uj)/N{Q) 
which is based on an MMP form with usf = 41 meV 
shows no distinct sharp structure. Thus, a relatively 
smooth spectrum extending over a large energy scale 
produces in the quasiparticle density of states only very 
small, gradual modulations which would be hard but not 
impossible to detect. 



IV. NUMERICAL RESULTS 

In Fig.^lwe show our numerical results for the real part 
of the quasiparticle self energy (only for lo > and 

made positive) as a function of frequency to based on nu- 
merical evaluation of Eqs. H16() for the renormalized Mat- 
subara frequencies. The black dotted curves apply to the 
normal and the black solid lines to the superconducting 
state. Both are at temperature T = 10 K. The spectrum 
Pxi^) used in the top frame of Fig.^jis the MMP form 
of Eq. with llIsf — 41 meV to which we have added a 
delta function contribution also a,t ue =41 meV. This is 
shown as the dashed line in the top frame of Fig. 13 The 
resulting boson structure in the normal state self energy 
is sharp and falls exactly at uje = 41 meV. It is the pres- 
ence of the 5-function in /^x(tj) which makes the peak 
so prominent. The corresponding peak in the supercon- 



ducting state (solid black curve) has shifted to higher fre- 
quency and falls slightly below to = uje + = 63.3 meV 
where Aq is the gap amplitude. For an s-wave supercon- 
ductor the shift in the boson structure would fall exactly 
at loe + Aq but for d-wave it falls below this value be- 
cause the gap is distributed in value and we are seeing 
the result of a distribution of shifts from to Aq. The 
peak has also broadened and the weight under it appears 
to have increased. In this sense, the opening up of the 
superconducting gap in the quasiparticle density of states 
has not only shifted the boson structure but has also, in 
a sense, enhanced it. 

This same statement applies even more strikingly to 
the corresponding optical quantities. For the imaginary 
part of the memory function, optical mass renormaliza- 
tion Aop(u;) = {[(rn*{uj)/m] — 1}, black dashed curve for 
the normal state and black dash-dotted curve for the su- 
perconducting state, the boson structure is much greater 
in the dash-dotted curve. Comparing the normal state 
memory function (black dashed curve) with the self en- 
ergy (black dotted curve) shows that the structure in the 
memory function is indeed at V2^£; rather than at uje 
although there is a background to our boson spectrum 
besides the (5-function. This square root of two shift be- 
comes critical in serious comparison ARPES and opti- 
cal results and has not been appreciated in the pastil 
Note that for the superconducting case the structure in 
the memory function is greatly enhanced over its nor- 
mal state value and that it is also shifted upwards in 
energy as for the self energy. Its exact position depends 
on the details of the d-wave gap structure and a complete 
Eliashberg calculation is required to determine it. The 
size of the peak in the underlying electron-boson spec- 
tral density I'^xi^) is not related simply to the size of 
the structure seen in the self energy or memory function. 

To produce sharp, easily identifiable, structures in 
ARPES or optical self energies, it is necessary to have 
correspondingly sharp structures in /^x(a;). This fact is 
well illustrated in the bottom frame of Fig. 0] where we 
show results based on a simple MMP model, Eq. H17|l . 
spectral density /^x(w) (solid curve in the top frame of 
Fig. O. Direct comparison of the curves in the bottom 
frame of Fig. ^ with the corresponding curves in the top 
frame shows that now the boson structures are much re- 
duced. Although there is a broad peak at 41 meV in 
/^x(a;) for the pure MMP spectrum, this does not trans- 
late into a discernable peak in any of our results in the 
bottom frame of Fig. ^ From this comparison we con- 
clude that a broad peak in Pxi^) is difficult to detect as 
a clear signature in either the real part of the quasipar- 
ticle self energy Si(a;) or in the optical mass renormal- 
ization wAop(a;) in the normal as well as in the supercon- 
ducting state. A detailed Eliashberg analysis is needed to 
extract J^x(w) from such data. A method for achieving 
this has been described in Refs. [lland[l3 An appropri- 
ate derivative of Tj^^(a;), namely (P[ivT~p^{u;)]/dLU^ which 
is closely related to Pxi^)^ see Eq. ©, is used to get a 
first modification to the MMP form, Eq. (|16|l . that fits 
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the optical data better. The procedure can be continued 
until a reasonable fit with the optical data is obtained 
from the model spectral density. However, no high accu- 
racy fit has so far been achieved. 

In the absence of such an analysis experimentalists 
have tried to get information on structure in Px{^) 
through comparison between normal and superconduct- 
ing state data. To this end we show in the bottom 
frame of Fig. 01 results for wAop(i^) in the normal state at 
T = 90K (black dash-double dotted curve). A subtrac- 
tion of this data from the dash-dotted curve in the super- 
conducting state gives the difference (gray dash-dotted 
curve). On comparing this curve with the input /^x(cl;) 
(solid line, top frame of Fig. |3Jl we see little detail cor- 
respondence. The peak in the difference curve is not at 
41 meV and the MMP tails at larger values of u> are not 
picked up. This holds as well if the normal state data at 
T = lOK (black dashed curve), albeit not accessible to 
experiment, is used in the subtraction (gray solid curve) 
or if Si(ijj) is considered (gray dotted curve). This sub- 
traction procedure does not give a reliable way to relate 
boson structure in experimentally measured quantities 
directly to boson structure in the electron-boson spectral 
density and is not recommended as an analysis proce- 
dure. This is also seen in the middle frame of Fig. ^ 
which applies to the case of the spectrum fit to optical 
data on YBCOe.gs twinned single crystals>i& It is shown 
as the dotted curve in the top frame of Fig. O It has 
a A = 2.01 and corresponds to an MMP form modified 
with an optical resonance peak at ojr = 41 meV and there 
is zero weight at small w. This spectrum is intermediate 
in "sharpness of boson structure" between the (5-function 
and the pure MMP case. We see that this reflected itself 
in how structured Si(a;) and Lo\op{uj) are. Note in partic- 
ular that the large peaks in the black dash-dotted curve 
for lo\op{lij) is now shifted to higher frequencies in the 
middle frame when compared to the top frame although 
the optical resonance peak is at the same frequency as 
is the (5-function in Fig. 13 (This is because the gap Ag 
is bigger now.) This demonstrates once more that no 
exact correspondence between size and position of the 
structure in bj\op{uj) and I'^xi'^) is possible without a 
detailed analysis, particularly for extended spectra. Also 
these functions strongly reflect the energy dependence 
of the underlying quasiparticle density of states that ex- 
ists in the superconducting state (dotted curve, bottom 
frame of Fig. 13) and it is not possible to, so to speak, sub- 
tract out the effects of these modulations from the data. 
Formally the subtraction procedure described above be- 
tween superconducting and normal state quantities when 
applied to the middle frame of Fig. ^ again fails to give 
reliable information on the shape of the boson spectrum. 

In the top frame of Fig.|Slwe show results for the optical 
scattering rate T~p{Lo) based on the model /^x(^) which 
consists of an MMP form plus a ^-function (dashed curve 
in the top frame of Fig. ^ at uje = 41 meV. The solid 
black curve is in the superconducting state at T = 10 K 
while the black dotted curve is in the normal state at the 



same temperature. While the underlying Pxi^) is the 
same in both cases yet the amount of structure seen in 
the solid curve is much more pronounced than it is for the 
dotted curve. This is due entirely to the energy depen- 
dence of the quasiparticle density of states (dashed line, 
bottom frame of Fig. ^ which arises in a d-wave super- 
conductor. On comparing the solid and the dotted curve 
we note a shift upward in the initial sharp rise in the scat- 
tering rate as we go to the superconducting case. This 
corresponds to the opening of the gap. The black dashed 
curve is in the normal state at T = 90 K. Temperature 
smears the sharp rise at small uj. Only the T = 90K 
data is available to experiment. If, as we did in Fig. ^ 
we take the difference between superconducting and nor- 
mal state optical scattering rate we get the gray solid line 
for T = lOK and the gray dashed curve for T = 90 K. 
The minimum in each of these curves is around 60 meV 
which is loe shifted by approximately Aq. It is obvi- 
ously hard to get reliable information on oj^; from such 
data. The middle and bottom frames of Fig. |S1 are for 
less structured I'^xi'-^) ^^'^ show a progressive decrease 
in the corresponding structures in r^^(w). 

We turn next to the exact correspondence between 
t~p(llj) and the quasiparticle scattering rates given as 
twice the imaginary part of the self energy. In the pre- 
vious section we have already examined this relation but 
based our discussion on approximations for the conduc- 
tivity formula. Next we use the more exact Eq. I|15|) and 
solutions of the full Eliashberg equations IjKil) . 



MODEL SELF ENERGY 



We return to the model function T^],j^f,i{i^) defined in 
Sec. mil It can be constructed from the numerical data 
shown in Fig. |Slfor T^p^{uj) defined as 



r-i(r,c^)-^Rea-i(r,^) 



and 



.[^) =^Ima-(r,c.) 
m J op 47r 



(Tl{T,Uj) 



47r af(T,Lo) + al{T,io)' 



(18) 



47r al{T,uj) + al{T,Lj) 



(19) 

We note in passing, that for small lo the optical effec- 
tive mass formula (|19() is dominated by the imaginary 
part of the optical conductivity <T2{T,lo) in the super- 
conducting state because it diverges as w"^ for small lo 
with coefficients proportional to the inverse square of the 
penetration depth. This is not related to a quasiparticle 
effective mass, but is rather a property of the condensate 
itself. Nevertheless, this is what has been done tradi- 
tionally and we need to follow this procedure here to 
make contact with the literature. An alternate approach 
would be to subtract out of the imaginary part of the 
optical conductivity the condensate contribution before 
forming the ratios indicated in Eqs. (|18|l and H19|) . In 
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this alternate approach the resulting ai and tT2 would 
refer directly only to the normal fluid part of the opti- 
cal conductivity and so would be more closely related to 
quasiparticle properties. But this is not what is done in 
the literature. 

For comparison with the quasiparticle self energy it is 
useful to use t~Ij^^i{T,lli) given by 



d 



Luai{T,uj) 



(20) 



than t^p^{T,ll>). Its imaginary part is given by the KK- 
transform 



a{Lu) 



1 



duj' 



(21) 



In Fig. we show results for t~1^^i{uj) based on the 
numerical results of Fig. O and compare directly with 
the imaginary part of the self energy T,2{uj) obtained di- 
rectly from Eliashberg calculations, namely Iml](w) = 
— Ima)(a; -I- iO"*") of Eq. ((TB|l . By definition of the self 
energy Cj{u! + iO~^) = uj — E(tj -f zO+). In the two top 
frames of Fig. we show results for the MMP form (solid 
line in the top frame of Fig. |3J) and in the two bottom 
frames the YBCOg.gs form (dotted curve in the top frame 
of Fig. |3J) for Px{'^) is used. In each case we present 
two frames, the left hand frame is for the normal state 
and the other frame for the superconducting state. All 
results are at temperature T = 10 K. For ease of com- 
parison with T~p{uj) (dashed hues) it is T~p{u))/2 (solid 

lines) and t^^j^^i{uj) / 2 (dotted lines) that is shown. Even 
when a factor of one half is included to make magnitude 
more comparable, the optical scattering rates differ sig- 
nificantly from the quasiparticle scattering rates. On the 
other hand, in the normal state the model scattering rate 
agrees almost perfectly with the quasiparticle scattering 
rate (top left hand frame). It is clear that t~1^i^i{uj) 
should be used in comparison with ARPES data and not 
T~p{uj). The close correspondence between model and 
optical rates is, however, lost when the superconducting 
state is considered (top right hand frame). The solid 
curve is significantly different from the dotted curve. In 
particular, T^],j^f,i{uj) is smaller at small w and rises faster 
around 50 meV after which it stays significantly above its 
quasiparticle counterpart up to almost 175 meV where 
the two curves cross again. These differences have their 
origin in the energy dependence of the superconducting 
quasiparticle density of states (solid line, bottom frame 
of Fig. 13) and have nothing to do with the boson struc- 
ture. For a comparison of Top^(w) (optics) and t^^(w) 
(ARPES) based completely on experiment, the reader is 
refered to Kaminski et alw^ 

When the underlying /^x(a;) is more structured, as 
is the case for the YBCOg.gs spectrum (dotted curve 
in Fig. 13), the correspondence between t~Ij^^i{uj) and 
t~p^{lu) is not as good as is seen in the two bottom 
frames of Fig. El There is a small additional oscilla- 
tion in the model case, not present in the self energy. 



Nevertheless, the agreement between them is still much 
closer than is the case for the optical rate (dashed curve) . 
Note, however, that for the superconducting state the dis- 
agreement between T~p^{uj) (solid curve) and t~1^^i{lj) 
(dotted curve) is now much greater and, in particular, 
a large peak is seen in the model curve which is not 
there in the dotted curve. This peak can be related to 
the broad shoulder in the quasiparticle density of states 
which follows the logarithmic singularity (dotted line, 
bottom frame of Fig. |3J). 



VI. CONCLUSION 

We have analyzed the relationship between the boson 
structure seen in the quasiparticle self energy (be it real 
or imaginary part) measured in ARPES experiments and 
the corresponding structure seen in the optical memory 
function determined in the infrared conductivity mea- 
surements. Starting first with the real part of the self 
energy a 5-function peak at w = w^; in the quasiparticle- 
boson spectral density I'^xi^) shows up as a logarithmic- 
like peak at UJ — oje in the normal state. In the supercon- 
ducting state a gap develops in the quasiparticle density 
of states and this introduces further structures in quasi- 
particle quantities. This effectively shifts to higher en- 
ergy the normal state boson structure, broadens it, and 
can make it appear more prominent. For an s-wave gap 
the shift would be Aq (gap amplitude) but for d-wave 
it is somewhat less because a distribution of gap values 
is involved. By contrast, for the memory function the 
boson structure in the optical mass renormalization is 
at \plujE in the normal state rather than at a;_E and it 
is much less prominent. Only a very small peak results 
which greatly reduces the value of such measurements for 
determining boson structure as was attempted recently 
in Ref . 0. Again, the boson structure shifts in the super- 
conducting state and can also appear more prominent 
as a result of the additional modulation brought about 
by the opening of a gap. These modifications, however, 
have nothing to do with the structures in t^xiuj). The 
two effects are not additive and require a full non linear 
Eliashberg analysis to disentangle in detail as we have 
provided here. When extended rather than (5-function 
spectra are used, the situation is even more complex. For 
example, a rather broad peak in ^ in the MMP 

form for spin fluctuations produces no peak at all, even 
in the superconducting state. 

So far we have described only the real part of the 
self energy and equivalent optical quantity. Much the 
same can be said about scattering rates. For a (5-function 
/^x(tj), the quasiparticle scattering rate jumps from zero 
to a finite value ai io = loe and remains unchanged af- 
ter that. On the other hand, the optical scattering rate 
starts from zero at to — uje and increases gradually to- 
wards the same finite value which it only attains at very 
high Lo. Thus, the signature of a 5-function is not singular 
as it is for the quasiparticle case. 
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We have also found that in the normal state a new 
model scattering rate can be introduced which we denote 
T~l,^^i{T,Lo) and define as d[u}T~j^{T,uj)]/dLL}. It is com- 
pletely determined from optical data and has the advan- 
tage that it follows much more closely the uj dependence 
of the quasiparticle scattering rate than does T~p(T,iLj) 
itself. Its KK-transform is very close to the real part 
of the quasiparticle self energy. We propose that it is 
this quantity that should be used in comparison between 
ARPES and optical experimental data. The comparison 
will be close in isotropic systems and in the anisotropic 
case, the model quantity gives an average over all direc- 
tions in the Brillouin zone. In the superconducting state 
no such simple comparison between ARPES and optics 
is possible. This needs to be kept in mind when analyz- 
ing experimental data, particularly in studies aimed at 
deriving boson structure from such data. This fact does 
not seem to have always been appreciated in analysis of 
experimental data, where features, at least partially as- 



sociated with the superconducting quasiparticle density 
of states, has been assigned to boson structure. 

The model spectral density Pxi^) consisting of an 
MMP form with superimposed the (5-peak was used to 
demonstrate that optics as well as ARPES just pick up 
the sharp structures and do not give direct information 
on a possible coupling to a smooth background. Never- 
theless, it may well be that it is precisely this coupling 
of the quasiparticles to such a background which is pri- 
marily responsible for superconductivity. 
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FIG. 4: Comparison of the real part of the quasiparticle self 
energy Ei(w) vs a; and optical effective mass renormalization 
LjAop(w) for the three model electron-boson spectral densities 
I^xi'^) shown in the top frame of Fig.|2l The top frame is for 
an MMP form with ujsf = 41 meV and \mmp ~ 0.554 aug- 
mented by a (5-function at ue = 41 meV with Xs = 0.946, the 
middle is obtained with a model I^xi^) obtained from con- 
sideration of optical properties of YBCOe.gs with A — 2.01 
(dotted curve in the top frame of Fig. and the bottom 
frame is for the pure MMP with lusf = 41 meV and A = 1.5. 
In all cases the curves come in pairs 'Ei{uj) in normal (black 
dotted curve) and superconducting state (black solid curve) 
and ujXop{uj) in normal (black dashed curve) and supercon- 
ducting state (black dash-dotted curve). In the middle frame 
we also include as (gray solid line) the difference between su- 
perconducting and normal state of ijj\op{ijj) and as gray dot- 
ted curve the difference between superconducting and normal 
state 'Si{uj). In the bottom frame similar difference curves 
are shown but with two different values of T in the normal 
state as labeled. 
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FIG. 5: Optical scattering rates T~p{T,uj) vs lo for the three 
quasiparticle-boson spectral densities used in Fig. 0] In each 
of the three frames the sohd black curve is for the supercon- 
ducting state at T = 10 K, the black dotted curve for the 
normal state at the same temperature, and the black dashed 
curve for the normal state at T = 90 K. The gray curves are 
the difference curves, solid curve for T = 10 K and dashed 
curve for T = 90K. 
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FIG. 6: Comparison of optical, Tgj^{u)) (dashed lines), quasi- 
particle, T~p^{uj) (solid lines), and model t^Ij^i{lo) (dotted 
lines) scattering rates vs lu for a temperature of T = 10 K. 
For easier comparison t~j^{ijj)/2 and /2 are shown. 

The two top frames compare our results for an MMP form 
with uJsF = 41 meV in the normal (left hand frame) and the 
superconducting state (right hand frame). The two bottom 
frames give the results for the model /■^x(<^) determined from 
optics for YBC06.95 as shown in the top frame of Fig. |3 (dot- 
ted curve). 



